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ABSTRACT 

The stability properties of newly born neutron stars, or proto-neutron stars 
(PNSs), are considered. We take into account dissipative processes, such as 
neutrino transport and viscosity, in the presence of a magnetic field. In order 
to find the regions of the star subject to different sorts of instability, we derive 
the general instability criteria and apply it to evolutionary models of PNSs. The 
influence of the magnetic field on instabilities is analyzed and the critical magnetic 
field stabilizing the star is obtained. In the light of our results, we estimate of 
the maximum poloidal magnetic field that might be present in young pulsars or 
magnet ars. 

Subject headings: convection - instabilities - MHD - stars: neutron - stars: 
magnetic fields 



1. Introduction 

Neutron stars are born in the aftermath of core collapse Supernovae explosions, the 
fate of massive stars within a range of 8 to SOM©. In the final stages of stellar evolution 
an iron core grows, until electron captures or photodesintegration of heavy nuclei trigger 
its gravitational collapse. The collapse proceeds until nuclear density is reached and the 
stiffening of the equation of state (EOS) provokes the bounce of the infalling material. A 
shock wave forms and propagates outwards, leaving behind a newly born, hot, lepton-rich 
PNS. These objects are prodigious emitters of neutrinos of all types, which dominate the 
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energetics of core collapse Supernovae and play a crucial role in the explosion mechanism 
by carrying away the binding energy from the PNS and depositing a part of it in the outer 
regions. 

Whether or not the energy relocated by neutrinos is able to relaunch the stalled shock 
is still a controversial issue (Burrows, et al. 2000; Rampp & Janka 2000; Janka 2000; Lieben- 
dorfer et al. 2001) that has led to speculation about alternative ways the explosion is powered 
such as convection, rotation, and magneto-hydrodynamics effects. In particular, it has been 
recognized by a number of authors that convective energy transport in the newly born PNS 
might play an important role by increasing the neutrino luminosity required to ignite the Su- 
pernova explosion (Burrows, & Fryxell 1993; Janka & Miiller 1996; Bruenn & Dineva 1996) 
Convection in PNSs can be driven not only by the lepton gradient, as originally suggested by 
Epstein (1979), but also by the development of negative entropy gradients, which is common 
in many simulations of supernovae models (Bruenn & Mezzacappa 1994; Bruenn, Mezza- 
cappa, & Dineva 1995; Rampp & Janka 2000; Liebendorfer et al. 2001) and evolutionary 
models of PNSs (Burrows & Lattimer 1986; Keil & Janka 1995; Keil, Janka & Miiller 1996; 
Pons et al. 1999, 2001a,b) despite differences in the equation of state and neutrino transport. 

Although stellar core collapse may also have dramatic consequences on the magnetic 
field strength, the effect of magnetic fields was not considered in the above mentioned works, 
mainly due to the tough technical difficulties inherent to multidimensional magnetohydrody- 
namical simulations. It is known that the magnetic field in a collapsing star can be amplified 
by many orders of magnitude due to the conservation of the magnetic fiux. This idea is often 
used to explain the origin of strong magnetic fields in neutron stars (Ginzburg 1964; Ginzburg 
& Ozernoy 1964; Woltjer 1964; Zeldovich & Novikov 1971; Shapiro & Teukolsky 1983). The 
conductivity of plasma inside the star is so large that the magnetic decay timescale is much 
shorter than the collapse time for a lengthscale of the magnetic field comparable to the stel- 
lar radius. Therefore, the field is frozen-in, and the magnetic flux is conserved during the 
collapse stage. In the case of a frozen-in magnetic field, the average field strength increases 
approximately as i? oc p^/^ where p is the density. By assuming that the magnetic field 
strength of the core of a massive star is comparable to that of strongly magnetized white 
dwarfs (~ 10^ G), at the end of core collapse, the neutron star may have a magnetic field 
~ 10^^ — 10^"^ G, comparable to that observed in young radio-pulsars. Following the same 
arguments, if we consider that the magnetic flux is conserved during the evolution of a mas- 
sive star from the main sequence, starting from an initial magnetic fleld of ~ 10^ — 10^ G, we 
end up with a magnetic field of ~ 10^^ — 10^^ G at the time the core of the star has become 
a neutron star. Similar arguments are often used to estimate the magnetic field of neutron 
stars formed from magnetized white dwarfs by the accretion induced collapse mechanism. 
Note that these estimates give the order of magnitude of a poloidal field. The toroidal field. 
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however, might be much stronger (Ardelyan, Bisnovatyi-Kogan, & Popov 1980). 

Under certain conditions, convective motions in PNSs can amplify the magnetic field 
via dynamo action. The most optimistic estimates lead some authors (Thompson & Duncan 
1993; Thompson & Murray 2001) to the conclusion that turbulent dynamo action could be 
responsible for the formation of ultra-magnetized neutron stars (magnetars), with fields as 
strong as ~ 10-*^^ G at a very early evolutionary stage. Certainly, such strong magnetic fields 
influence magneto-hydrodynamic processes in PNSs, particularly their stability properties. 
At this point, the fact that the convective overturn is intrinsically a 3D phenomenon makes 
fully consistent numerical simulations a very laborious and computationally expensive way 
to address the problem, and simplified semi-analytical approaches are extremely useful to 
evaluate the necessity of more detailed studies. In this respect, some attempts have been 
made to establish definite criteria as to the different types of instability and their growth 
times (Bruenn & Dineva 1996; Miralles, Pons, & Urpin 2000, MPU in the following). In 
MPU we derived the general criteria for instability, taking dissipative processes into account 
such as neutrino transport and viscosity. In the present paper, we consider the stability 
properties of magnetic PNSs in detail. It is well known (see, e.g., Chandrasekhar 1961) that 
the magnetic field may stabilize a fluid against convection causing an additional effective 
viscosity due to Lorentz force. On the other hand, the field can lead to new branches of 
instability associated, for instance, with the Alfven waves. In this paper, we present the 
main results of a general stabihty analysis in magnetized PNSs. The stability properties will 
be considered for a wide range of magnetic fields. 

The paper is organized as follows. In §2, we obtain the dispersion equation which deter- 
mines the condition of linear instability from the equations governing magneto-hydrodynamics 
of PNSs and neutrino transport. In §3, we discuss the stabihty properties of both non- 
magnetic (§3.1) and magnetic (§3.2) PNSs, apphng the criteria of hydromagnetic stability 
to the results of numerical simulations to identify the unstable zones. We also calculate the 
critical magnetic field that stabilizes a PNS. Finally, in §4 we summarize and discuss our 
main findings. 

2. The dispersion equation 

Consider the stability properties on adopting a plane-parallel geometry. We assume 
that within the layer between z — and z — d the gravity g is directed in the negative 
z-direction and neglect a non-uniformity of g as well as general relativistic corrections to the 
equations of hydrodynamics. The characteristic cooling time scale of a PNS is assumed to 
be much longer than the growth time of instability thus the latter can be treated in a quasi- 
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stationary approximation. Since the convective velocities are typically much smaller than 
the speed of sound, one can describe instability by making use of the standard Boussinesq 
approximation (see, e.g.. Landau & Lifshitz 1987). We consider the linear instability regime 
when the equations governing small perturbations can be obtained by the linearization of 
hydrodynamic equations. In the following, small perturbations of hydrodynamic quantities 
will be marked by a subscript "1". For simplicity, the unperturbed magnetic field is assumed 
to be vertical. The linearized momentum, continuity and induction equations read 



pvi = -Vpi + gpi + —(V X Bi) X B + pi/Avi, (1) 

V • vi = 0, (2) 
Bi = V X (vi X B) + i/„ABi, (3) 

V • Bi = 0, (4) 



where p and p are the pressure and density, respectively, v is the viscosity and is the 
magnetic diffusivity. We assume that the matter inside a PNS is in chemical equilibrium, 
thus the density can generally be considered as a function of the pressure p, temperature 
T and lepton fraction Y, Y — {rie + n^)/n, with rie and rii, being the net (particles minus 
antiparticles) number densities of electrons and neutrinos, respectively, and n — Up+Un is the 
number density of baryons. Then, the perturbations of density, entropy per baryon (s) and 
neutrino chemical potential (fi) can be expressed in terms of pi, Ti and Yi. In the Boussinesq 
approximation, the perturbations of pressure are negligible because the fluid motions are 
assumed to be slow and the moving fluid elements are nearly in pressure equilibrium with 
their surroundings. Therefore, we have 

Pi ^ -p + SY,^ , (5) 

si ^ mpCpj; + aYi, (6) 

where rUp is the proton mass (we neglect the mass difference between protons and neu- 
trons), f3 and 6 are the coefficient of thermal and chemical expansion, respectively, f3 = 
-{d\np/d\nT)pY, 6 = -{d\n p/dY)pT; a = {ds/dY)pT; and Cp = {T/mp){ds/dT)pY is the 
specific heat at constant pressure. 

The above equations should be complemented by the equation driving the evolution of 
chemical composition and heat balance. We employ the equilibrium diffusion approximation, 
which is sufficiently accurate and reliable during the early stage of PNS evolution, when the 
mean free path of neutrino is short compared to the density and temperature length scales. 
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In this approximation, the diffusion and thermal balance equations can be linearized (see 
MPU for details) and written as 

AT" 

+ vi • Vr = + AyAFi, (7) 

fi AVT ATi 

^ - Vi • = + KyAY,, (8) 

where 

AVT = (Vs-aVY)={ — \ Vp-VT (9) 

mpCp \op J s,Y 

is the superadiabatic temperature gradient, i.e., the difference between the temperature gra- 
dient of a fluid with constant entropy and composition and the actual temperature gradient, 
and the transport coefficients Xt, kt, Ay, /ty have been defined in MPU. 

Equations (l)-(4), (7) and (8) together with the corresponding boundary conditions 
determine the behaviour of small perturbations. For simplicity, we consider the case of 
vanishing perturbations at the boundaries z = and z = d. Note that other boundary 
conditions cannot change the main conclusions of our analysis qualitatively. 

Without loss of generality, the dependence of all perturbations on time and on the 
horizontal coordinate can be chosen as exp(7t — ikx), where k is the horizontal wavevector 
and 7 is a constant which can be complex and its real part gives the inverse growth (or decay) 
timescale. Then, equations (l)-(4), (12), and (13) can be reduced to only one equation of 
higher order. Solving for viz, we obtain 

A[(7 - «;tA)(7 - Ay A) - A^KyA^] 



(7-^A)(7-i/^A)-4^^2 



Viz 



AWT 



( dY 

[/?(7 - Ay A) + 5AtA] I (7 - VmA)viz, (10) 

where = B / ^JA-Kp is the Alfven velocity and A = ^ — k'^- The coefficients of this equation 
are constant in our simplified approach. The solution for the fundamental mode with "zero 
boundary conditions" has the form viz oc sm{nz/d). Thus the dispersion equation for the 
fundamental mode is 

[(7 + (jJ^) (7 + UJrri) + "^l] [(7 + ^t) (7 + ^y) - LOtY^yt] = 
il + ^m) ioli'l + i^Y - OIUyt) + -^^^^ , (H) 
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where a — 8/ (5. In this expression, we introduced the characteristic frequencies 

ujyt = ^tQ"^, ooty = hyQ^, ooa = CAir/a, (12) 

~ TQ^ ' '^■^ ~ ~~Q^ ■ d^' 

with = (vr/rf)^ + k"^. The quantities ou,^, ou^, out, and uy are the inverse time scales of 
dissipation of the perturbations due to viscosity, magnetic diffusivity, thermal conductivity 
and chemical diffusivity, respectively; Uyt characterizes the rate of diffusion caused by the 
temperature inhomogeneity (thermodiffusion), and ujty describes the influence of chemical 
inhomogeneities on the rate of thermal evolution; ujg is the frequency (or, in the case of 
instability, the inverse growth time) of the buoyant wave; ul characterizes the dynamical time 
scale of the processes associated with the lepton gradient, and cua is the Alfven frequency. 



3. Stability of proto- neutron stars 

3.1. Non-magnetic PNSs 

Consider initially the case when B = 0. Then, equation (11) is reduced to a cubic 
equation with the roots given by 

7^ + a27^ + ai7 + ao = 0, (13) 

where 

0,2 = COu +(^T + ^Y, 

tti = UtOJy — ^TY^YT -\- ^v{^T -\- I^y) — {l^l -\- 
Oq — iiJ uitiJj'ijJY — Uxy^Yt) 

— oj'1{ujy — (xujyt) — <^L i'^T <^ty] ■ (14) 



a 



Equation (13) describes three essentially different modes, corresponding to the three 
roots of the equation, which can generally exist in a chemically inhomogeneous fluid. The 
condition that at least one of the roots has a positive real part (unstable mode) is equivalent 
to one of the following inequalities 



a2 < 0, ao < 0, aia2 < ao 



(15) 
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being fulfilled (see, e.g., Aleksandrov, Kolmogorov, & Laurentiev 1985; DiStefano III, Stub- 
berud, & Williams 1994). Since kt and Ay are positive defined quantities, the first condition 
a2 < will never apply, and only the two other conditions determine the stability of non- 
magnetic PNSs. 

Generally, the last two conditions (15) are rather complex and depend on the horizontal 

wavevector of perturbations, k. The critical temperature and lepton gradients which deter- 
mine the transition from stability to instability can be quite different for perturbations with 
different k. In order to obtain the true criteria wc have to find the minimal values of these 
gradients which can cause instability. From the properties of the kinetic coefficients we have 
{ktXy — kyXt) > 0, therefore the gradients in the last two conditions (15) are minimal for 
perturbations with the wavevector k minimizing the quantity Q^/k'^. This happens when 
k"^ — {'K/dy /2. For such a horizontal wavevector, the instability criteria read 

^{5kt - (3ky) - {(3Xy - SXt) 
az 1 

27 TT^ 

< ~ '^y^t), (16) 

dY AVT 

-^[5{iy + \y) + Pky] - + kt) + 5\t] 

< -^(«r + Ay)[(i/ + KT)(i^ + Ay) - ArKy]. (17) 

Conditions (16)-(17) divide the AVT — Vy plane into four regions, which are charac- 
terized by different stability properties. The size and configuration of these regions depend 
on the thermodynamic and kinetic properties of nuclear matter. The detailed analysis of 
stability of non- magnetic PNSs based on criteria (16) and (17) has been done in MPU. In 
that paper, as in the present one, we use results from numerical simulations of PNS evolution 
by Pons et al. (1999) to calculate the different thermodynamical derivatives, diffusion coeffi- 
cients, and conductivities appearing in the stability criteria. To illustrate the large variety of 
situations encountered during the Kelvin-Helmholtz phase of a PNS, we have chosen three 
different cases. The physical conditions and the different transport parameters are summa- 
rized in Table 1. Case A corresponds to a layer near the center of the star, with an enclosed 
mass of 0.05 Mq, at an evolutionary time of t = 2 s, when the PNS is in the deleptonization 
stage and the core is being heated. Cases B and C correspond to two different layers at 
a time of t = 20 s, during the cooling stage, mainly driven by thermal neutrinos since an 
important part of the lepton content of the star has already been radiated away. The layer 
corresponding to case B is located at the same depth as A (0.05 M©) star while the layer in 
case C has an enclosed baryonic mass of 1 M©. In our calculations we have used the value 
of d given by the pressure height d — l= \ d In p/dr \~'^. 
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In Figure 1 we plot the lines of critical stability given by equations (16) and (17) (solid 
and dotted line respectively) in a AVT — plane. The asterisks indicate the particular 
values of the gradients for each of the three cases considered, the parameters of which are 
summarized in Table 1. The shaded area covers the region where all the roots of the cubic 
equation (13) are real. We denote by neutron fingers a region where condition (16) is fulfilled 
but not condition (17); by convection we denote those regions where both conditions (16) 
and (17) are satisfied; we denote by semiconvection the case when condition (17) is satisfied 
but condition (16) is not. 

By applying Routh critcrium (DiStefano 111, Stubbcrud, & Williams 1994), which states 
that the number of unstable modes of a cubic equation is given by the number of changes of 
sign in the sequence 



we deduce that the number of unstable modes in the regions labelled by convection or neutron 
finger is always one. Thus, the root corresponding to this mode is real and positive. Notice 
that, while in the convection case all roots are real, in the neutron finger case, the two stable 
modes can be either real or complex conjugates. In the region labelled semiconvection we 
have two unstable modes, which may be either real or complex conjugates. Semiconvection 
can be, therefore, oscillatory (complex conjugate unstable modes) or non-oscillatory (two real 
unstable modes). We emphasize that the convention we use to divide the plane AVT — VY 
according to the instability criteria (16) and (17) differs from that of Bruenn & Dineva (1996). 
According to their convention, the convectively unstable region is extended to contain that 
part of our neutron finger and semiconvection regions where all roots arc real. In general, the 
classification in different types of instability is somehow arbitrary, because there is no clear 
qualitative distinction between neutron fingers and convection: in both cases the unstable 
root is real. Our classification simply responds to the fact that, in the limit of vanishing 
diffusivity, criteria (16) and (17) reduce to the classical Rayleigh- Taylor and Schwarzschild 
criteria, and certainly is not the only valid choice. The important idea to bear in mind 
is that, when all the kinetic coefficients are of similar magnitude, the distinction among 
the different types of instabilities is less evident and only the concept of overall stabihty is 
meaningful. 




(18) 



3.2. Magnetic PNSs 



In the general case, B 7^ 0, equation (18) can be rewritten as 



7^ + ?37^ + ?27^ + ?i7 + ?o = 0, 



(19) 
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where 

gs = 02 + CUm, 

?2 = ai+ L0ma2 + 

qi — ao + cUmtti + u;\{u;t + u;y) 

go = cumao + u}\{utouy - outyooyt) ■ (20) 

Equation (19) describes four essentially different modes which can generally exist in chemi- 
cally inhomogeneous fluid in the presence of the magnetic field. The Routh criterium applied 
to a quartic equation states that the number of roots with positive real part (unstable modes) 
is given by the number of changes of sign in the sequence 

Jt m2-qi g|go 1 

I g3 qm - qi J 

Since ^3 is always positive, there is, at least, one stable mode. The criteria for hydro- 
magnetic instability is equivalent to requiring that any one of the following inequalities be 
satisfied: 

, (22) 
, (23) 

. (24) 

In general, the magnetic field has a stabilizing effect on all three conditions. Consequently, 
there will be a critical value of B above which none of the instability conditions is fulfilled 
and all modes become stable. 

To warm up, we begin analyzing the behaviour of the roots with increasing magnetic 
field. In Figure 2 we draw a sketch, plotting the real part (7) of the four modes as a function 
of u}\. Modes with positive 7 are unstable. In order to show more clearly the qualitative 
effect, we have assigned values to the dissipative frequencies which do not correspond to real 
values encountered in the early stages of a neutron star's life. The qualitative behaviour 
using real coefficients is very similar, although the different scales would not allow one to 
capture the trends in a simple plot. In the limit of a vanishing magnetic field {lva 0) 
there are three non-magnetic modes (sohd lines), discussed in the previous section, which 
remain essentially unaltered, and a stable magnetic mode with a value equal to —cUm (dashed 
line). The situation shown in Figure 2 would correspond to the convection case (1 unstable 
and 2 stable non-magnetic modes when B = 0). By switching on the magnetic field, the 
real part of the new magnetic mode (dashes) increases, and becomes positive in the point 



(21) 



go < 



gi 



g3g2 - gi < 
gjgo ^ 
(g3g2 - gi) 
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labelled a, which corresponds to the physical conditions such that qo — 0. At that point, the 
original unstable mode (upper solid line) has been partially damped. Further increasing the 
magnetic field, the real parts of the originally unstable mode and the new unstable magnetic 
mode merge, to form a pair of complex conjugate modes (point b). Convection becomes, 
thus, oscillatory, due to the stabilizing action of the magnetic field. At a critical magnetic 
field, given by condition (24), the real part of the unstable modes vanishes (point c), and 
for higher magnetic field all modes are stable. The originally stable modes are not seriously 
affected. 

The general trends discussed above (namely: damping of the unstable mode, a new 
magnetic mode becoming unstable, merging of modes and change to oscillatory convection, 
and final suppression of instabilities) are also found when realistic kinetic coefficients from 
numerical simulations are used. To quantify the characteristic values of the magnetic field 
causing the effects described above, we have taken the parameters corresponding to case B 
from Table 1. Those modes with a typical wavelength of the order of the pressure height (see 
below for further discussion), suffer the transitions a, b, c in Figure 2 when B = 2.5 x 10^^ G, 
1.8 X 10^^ G and 2.09 x 10^^ G, respectively Thus, a field larger than 2.09 x 10^^ G entirely 
suppresses all types of instability in the considered region. 

In order to understand the effect of the magnetic field on each instability condition, it 
is more convenient to rewrite equations (22-24) separating the terms containing the Alfven 
frequency and magnetic dissipation: 



These inequalities can be further simplified if we take into account that the frequencies 
entering these conditions are of different orders of magnitude in PNSs. The dissipative 
frequencies, u^, Ut, ujy, ujty, and Uyt are approximately comparable to each other and 
range between 10 — 10^ s~^. The magnetic dissipative frequency, uOm, is much smaller. 
The electrical resistivity of hot nuclear matter (7^) is relatively low, TZ ^ 6 x 10~^^T8^ s 
(Yakovlev & Shalybkov 1991), thus the characteristic time scale of dissipation of the magnetic 
field is extremely long compared to other dissipative time scales. Taking into account that 
i^m — c^7?./47r, we can estimate LUm ~ 10~^^ s"^ for k ~ n/i. Therefore, we can neglect 
cUm when it enters conditions (23) and (24) in a row with other dissipative frequencies. In 




(25) 



(26) 



[ao + uJmai + oj\{ujt + a;y)][a2ai - ao + 0Jma2{(i2 + ojm) + a;l(<^,. + a;^)] 

< (a2 + UJmf X [{'^mO'O + Uj\{uJt<^Y - UtY<^Yt)\- 



(27) 
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addition, due to the small magnetic diffusivity, it also turns out that cum^^i, <^ cu^ and 
equation (22) can be written as 



—cuKcuy — CKCUyt) — ^ti^T <^ty \ H —{utcoy — uty<^yt) < 0. (28) 

All three conditions depend on the wavevector k. It is well known (Chandrasekhar 1961) 
that in the inviscid limit the minimal values of gradients which can cause instability of a 
magnetized fluid correspond to A; — > oo. Taking account of viscosity increases the most 
unstable horizontal wavelength and the critical gradients are reached for a large, though 
finite, value of k. In the case of low magnetic dissipation, k is much larger than ir/d. 
Substituting k ^ Q into equation (35), we have 

^{Skt - Pky) - ^^{PXy - SXt) < --r^ — (i^tXy - kyXt). (29) 

az 1 Agpa'^Urn 

Due to the particular physical conditions in PNSs {uOm ^ 1), even a relatively low magnetic 
field can make this condition not satisfied. We must point out, however, that this criterium 
is associated to the change of sign of new the magnetic mode (dashed line in Figure 2), which 
is irrelevant, because the dominant unstable mode is always that associated to one of the 
other two conditions. 

Let us turn back to the other two conditions, (26) and (27). The dissipative frequencies 
are usually small compared to the dynamical frequencies, ujg and wl, for typical length- 
scales in a PNS. For example, the contribution of terms proportional to the cube of dis- 
sipative frequencies in equations (26)-(27) does not exceed 10%. Therefore, we can ne- 
glect the terms containing only dissipative frequencies compared to those containing the 
dynamical ones (including the Alfven frequency). This latter approximation is justified 
if the magnetic field is large enough for uja to be larger than any dissipative frequency, 
uja > cudiss = Taax{uj,y,ujT,ujY,^TY,^YT)- Siucc cjdiss ~ 0.1 — 1 s~^, the condition uja > t^diss 
is fulfilled if the magnetic field exceeds some characteristic value which is of the order of 
10^^ — 10^'^ G. For B < 10^'^ G the infiuence of the magnetic field is insignificant and condi- 
tions (26)-(27) are exactly equivalent to conditions (16)-(17). For B > 10^^ G, the instability 
conditions can be simplified as follows: 

—UJg{uJY — OIIjJyt) — {^T — ^ + '^K'^T + <^y) < 0, (30) 

—uj'^g{ujT + 0Jv + auYr) - {^y + + + ^u^\ < 0. (31) 

As for equation (28), the minimal values of the gradients which cause the instability corre- 
spond to k ^ Q (short wavelengths). Substituting k — Q into equations (30) and (31), we 
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have 

^(5kt - (3ky) - ^(/?Ay - 5Xt) < -^^{i^t + Ay), (32) 

^[S{U + Xy) + P^y] - ^[Piu + Kt) + 5Xt] < (33) 

Note that these criteria are analogous to the instabihty criteria for non-magnetic PNSs (16)- 
(17), but replacing 

u{KTXy - KyXr) ^ ^^^3^ (^T + Ay), (34) 

(PB^ 

{kt + Xy) [{ly + Kt) {ly + Xy) - XrKy] ^li^'^p ^' ^^^^ 

and the magnetic field plays the role of an efficient dissipative process. In magnetic PNSs 
the dissipative effect caused by the Lorentz force is much greater than neutrino transport or 
neutrino viscosity effects. 

In Figure 3 we show the unstable zones at each time of the evolution of a young PNS from 
one of the numerical simulations in Pons et al. (1999). Top and bottom panels correspond 
to criteria (32) and (33), respectively. Criterion (33) depends on the neutrino viscosity, for 
which available estimates are sometimes contradictory (van den Horn & van Weert 1981; 
Goodwin & Pethick 1982; Thompson & Duncan 1993). In our simulations we used the 
value V — rfKT (van den Horn & van Weert 1981). The contours indicate the convectively 
unstable region for the following values of a passive magnetic field: 10^^ G (sofid), 5 x 10^^ 
G (dotted), 10^6 G (dashed), 3x 10^^ G (dash-dot) and 5 x 10^^ (dash-dot-dot-dot). It turns 
out that, with this viscosity, the critical magnetic fields for criteria (32) and (33) arc of the 
same order of magnitude. The magnetic correction on (33) is somewhat more restrictive 
and, generally, its associated instability is suppressed with weaker magnetic fields. The 
difference is caused mainly by the different combinations of kinetic coefficients that appear 
in the equations. Both criteria require a rather strong magnetic field, S ~ (1 — 5) x lO^^G, 
in order to stabifize the fluid completely, although the region where condition (33) is fulflUed 
is typically smaller than that given by condition (32). Note that moderate magnetic flelds 
in the order of < 10^^ G leave the unstable regions essentially unchanged, and the effects 
are visible only for B > 10^^ G. 

In Figure 4 we plot the critical magnetic held (-Bgt) which completely stabifizes a certain 
layer of the PNS (corresponding to enclosed baryonic masses of 0.2, 0.5, 1.0, and 1.4 Mq) 
as a function of time. Layers with Mb < IMq are stable at the beginning of the evolution. 
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Later on, when the corresponding layer becomes unstable, Bst grows rapidly, reaching some 
maximum value, and then decreases monotonously as the temperature and lepton gradients 
are smeared out via neutrino diffusion. The maximum value of B^i at different points varies 
in a very narrow range (2 — 5 x 10^^ G), but this maximum is reached at different instants in 
time. Note that -Bst provides stabilization in the largest unstable lengthscale. Stabilization 
at smaller scales is accomplished with weaker magnetic fields. The instability seems to be 
especially persistent in the central region with Mb < 0.5Mq, where convective motions might 
last for a long time (Keil, Janka & Miiller 1996) and ~ 5 x 10^^. 



4. Discussion 

In this paper wc have considered the effect of magnetic field (consistently coupled to 
other dissipative processes, such as neutrino transport and neutrino viscosity) on the convec- 
tive instability of PNSs. Our results show that transport phenomena play an important role 
in determining the size of the region of the star which is subject to instability. The particular 
values of the transport and kinetic coefficients also determine the kind of instabilities which 
can be present in an unstable region {convection, neutron fingers, semiconvection) . 

For low magnetic fields {B < 10^^ G) the region where instabilities are present remains 
unaffected when compared to the non-magnetic case. That region covers an important part 
of the PNS and lasts for almost the whole Kelvin-Helmholtz phase of the PNS evolution. 
Gonsequently, due to the existence of such an extended region subject to instabilities, con- 
vective motions are likely to appear and turbulent dynamo action can significantly increase, 
locally, the magnetic field strength. 

Once the magnetic field has a significant strength (10^^ G < B < 10^^ G), its effect 
is twofold. First, it is responsible for the appearance of a new mode which can generally 
be unstable, although in this case, its growth time is always longer than that of the other 
unstable modes. Secondly, the magnetic field acts as a stabilizing agent for all types of 
instability studied in the non-magnetic case, progressively reducing the size of the convective 
regions. Our calculations show, however, that only a relatively strong magnetic field can 
infiuence significantly the stability properties of PNSs, particularly in the central region with 
Mb < 0.5Mq. For instance, if the dynamo action is sufficiently efficient, the duration of the 
convective epoch can be noticeably shorter as the magnetic field increases up to > lO^^G. 
We find that the critical magnetic field that suppresses all sorts of instability happens to be 
in a narrow range, at about 5 x 10^^ G. This can be considered as an upper limit to the 
poloidal magnetic field generated by convective dynamo. 
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In our simplified approach, we have considered the case of a vertical magnetic field. 
It is known (Chandrasekhar 1961) that this magnetic configuration requires the maximal 
field strength for stabilization. In real conditions, the field geometry may be more complex 
with a substantial component perpendicular to the temperature and composition gradients. 
In PNSs, this field component can correspond, for example, to a toroidal magnetic field. 
For this geometry, the stabilizing magnetic field is typically several times lower. Therefore, 
stabilization against hydromagnetic instabilities can really take place at a slightly weaker 
magnetic field than that obtained in our calculations. In the presence of differential rotation 
(a very likely situation after the core collapse), the main contribution to the total field 
strength would be given by the toroidal magnetic field (St) and the poloidal component 
{Bp) would be less important. In a quasi-stationary regime, the rate of stretching of the 
toroidal field lines from the poloidal configuration by differential rotation is compensated 
by turbulent dissipation (Parker 1979). By equating both rates one can write the following 
equation: 

vriT^,^^^ (36) 



r 3 

where r is the characteristic stretching time, vt the turbulent velocity and £t the turbulent 
lengthscale. Since the stretching is due to differential rotation, r might be estimated as 
r =1 AQ 1^^, AQ being the difference between the maximum and minimum angular velocity. 
An estimation of V^-Bt is ~ Bt/ R^, R being the typical radius of the star. Using the following 
values, ~ 20 km, vt ~ 10^ cm/s and £t ~ lO^cm, we obtain 



Bt AQ 

Wr, ^ li^ 



(37) 



p 

,3 o-l^ 



Assuming that neutron stars are born with large angular velocities (say Q = 10 s~ ), and 
that AQ/Q is comparable to that of the Sun (0.1 — 0.3), we can estimate the toroidal field to 
be 100 — 300 times stronger than the poloidal component. Since the total magnitude of the 

magnetic field generated by dynamo action should not exceed, in any case, 10^^ G we have 
to conclude that the poloidal component which may represent the observed pulsar field must 
be in the range 3 x 10^^ — 10^"^ G. This estimate is in agreement with the observed fields of 
young pulsars in supernova remnants. The toroidal field could be, however, much stronger, 
which might lead to important consequences on the structure of neutron stars (Bocquet et 
al. 1995; Cardall, Prakash, & Lattimer 2001). 
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Table 1: Physical conditions and transport coefficients for three representative situations 
during the first seconds of life of a PNS. Cases A,B,C represent the three possibihties found 
in numerical simulations, classified according to criteria (16)-(17). 



Case 


Thermodynamic state 


Transport coefficients (cm) 




ub (fm-3) T (MeV) Yl (3 5 


Kj' Ky Aj" Ay 


A (stable) 

B (convection) 

C (fingers) 


0.301 21.1 0.333 0.103 0.111 
0.350 33.6 0.259 0.145 -0.058 
0.197 12.1 0.096 0.069 -0.570 


0.556 -0.851 -0.046 0.496 
0.436 -0.536 -0.038 0.285 
2.46 -20.1 -0.23 4.67 
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Fig. 1.- 

Stability plane for the three thermodynamic states detailed in Table 1. Four different 
regions in the plane are separated by criteria (16)-(17) (solid and dotted hnes, respectively). 
The shaded area indicates the region where the three roots of equation (13) are real. The 
particular values of the gradients in each studied case (A,B,C) are indicated by asterisks. 
The scale factor I has been taken as the pressure height of each case. 

Fig. 2.- 

Sketch showing the behaviour of the real part of the roots of equation (19) as a function 
of the square of the Alfvcn frequency, for fixed values of the dynamical and dissipative fre- 
quencies. Dashed line corresponds to the mode decoupled in the non-magnetic case (referred 
to as magnetic mode in the text). Units are arbitrary and the points a, h and c indicate 
changes of sign and character (real-complex) of the modes. See §3.1 for further details. 

Fig. 3.- 

Contour plots showing the variation of the different convectively unstable zones, assum- 
ing different strengths of the magnetic field. The contour lines correspond to the following 
values of the magnetic field: 10^^ G (solid), 5 x 10^^ G (dotted), 10^^ G (dashed), 3 x 10^^ 
G (dash-dot) and 5 x 10^^ (dash-dot-dot-dot). The top panel refers to the unstable region 
according to criteria (32) and the bottom panel refers to (33). The x-axis indicates the 
evolution time and the y-axis the enclosed baryonic mass of the layer. The innermost part 
of the star remains unstable at intermediate times even when magnetic fields as large as 10^^ 
G are present. The asterisks correspond to the cases A, B and C described in the text. This 
PNS model was taken from numerical simulations performed in Pons et al. (1999). 

Fig. 4.— 

Critical magnetic field strength that completely stabilizes the fiuid as a function of the 
evolution time. The four lines stand for four different layers, with enclosed baryonic masses 
of 0.2 (sohd), 0.5 (dotted), 1.0 (dashed) and 1.4 (dash-dot) M©, respectively The PNS 
model is the same as in Figure 3. In general, the deeper is the layer, the larger magnetic 
field is required for unconditional stabilization. 
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